Spin-charge separation in a strongly correlated spin-polarized chain 
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We combine the first-quantized path-integral formalism and bosonization to develop a phenomeno- 
logical theory for spin-charge coupled dynamics in one-dimensional (ID) ferromagnetic systems with 
strong interparticle repulsion, at low temperatures. We assume an effective spin-charge separation 
and retain the standard Luttinger-liquid plasmon branch, which is explicitly coupled to a ferromag- 
netic spin- wave texture with a quadratic dispersion. The dynamic spin structure severely suppresses 
the plasmon peak in the single-particle propagator, in both fermionic and bosonic systems. Our 
analysis provides an effective theory for the new universality class of ID ferromagnetic systems, 
capturing both the trapped spin and propagating spin-wave regimes of the long-time behavior. 
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Introduction. The Tomanaga-Luttinger liquid (TLL) 
is one of the most established examples of a non-Fermi 
liquid in low dimensions. The TLL theoretical frame- 
work captures many essential features of a large class of 
correlated one-dimensional (ID) quantum systems,— both 
fermionio^ and bosonic^ Recent theoretical investiga- 
tions have also focused on interacting ID systems that lie 
outside the conventional TLL paradigm^&£& & 10 ' 11 i 12 
At lower densities when the Coulomb repulsion dom- 
inates, the system becomes reminiscent of a fluctuat- 
ing Wigner solid.— A regime that has recently attracted 
a considerable interest is when the spin exchange cou- 
pling J is exponentially suppressed by strong repulsions, 
J <C T <C Ep, where Ep is the Fermi energy4&i2 
Explorations into this spin-incoherent limit began with 
the pioneering Bethe ansatz calculation of Cheianov 
and Zvonarev^^ who showed that the single-particle 
Green's function contains exponential decay in the spin 
sector instead of the usual TLL power-law singulari- 
ties in both spin and charge sectors. This makes this 
regime even more universal than the standard TLL, 
since the spin sector is fully thermalized. Moreover, the 
power-law exponents admit a nonunitary conformal and 
non-TLL classification^ Additionally, others have also 
looked for non-TLL behavior in ID systems by explor- 
ing the band curvature effects&i£i beyond the custom- 
ary linearized dispersions, as well as finite-energy effects 
stemming from an exact treatment of the spin sector, 
not relying on the bosonization of spin modes Very 
recently,— Zvonarev et al. explored low-energy prop- 
erties of strongly-correlated spin-polarized Bose gases, 
pointing out the existence of two regimes for spin prop- 
agation: a trapped spin regime, when spin flips are un- 
likely and spin excitation is mainly carried by charge fluc- 
tuations, and the open (propagating) spin-wave regime, 
where the spin excitations are carried by ferromagnetic 
magnons. Here, we develop an effective theory which 
helps to clarify the underlying physics and to generalize 
the results of Ref. Il3l 

In this paper, we continue in the spirit of the earlier 
investigations and calculate the single-particle Green's 
function G(x, t) of a ID spin-polarized chain of either 



bosons or fermions, focusing on the quadratic spin- wave 
curvature effects. We compare the different asymptotic 
forms of G(x, t) for both finite and vanishing exchange 
coupling J. The standard TLL model is antiferromag- 
netic, in accordance with the Lieb-Mattis theorem, 14 
which forbids spontaneous magnetization in a purely ID 
quantum system. Therefore, a spin-polarized correlated 
system is inherently outside the scope of the conventional 
TLL theory. Others have discussed a possible circumven- 
tion of the Lieb-Mattis theorem in ID systems with weak 
itinerant magnetism by adding a finite width to make the 
system effectively quasi-lDi 15 ' 16 These authors have also 
proposed that a quasi- ID spin-polarized chain is the can- 
didate system responsible for the (0.5 — 0.7) x 2e 2 /h con- 
ductance anomaly. There are also strictly ID fermionic 
models manifesting itinerant ferromagnetism that fall 
outside the class of systems considered by the Lieb-Mattis 
theorem^ One could also envision optically trapped ul- 
tracold ID Bose gases 1 ^ that become ferromagnetic in 
the presence of strong interparticle repulsion. 

In our system, the band curvature effects are present 
in the quadratic dispersion of the ferromagnetic spin- 
wave excitations. We will assume a strongly interacting 
Wigner crystal-like limit with individual particles having 
very low probability of exchanging positions with their 
neighbors. 7 The residual probability of tunneling through 
the interparticle repulsive potential gives rise to a small 
exchange coupling, which characterizes low-energy spin 
excitations. In the spin-incoherent TLLj 1 ^ the spin ex- 
citations are assumed to be fully thermalized by a fi- 
nite temperature. Here, we will take the temperature to 
be below the characteristic exchange coupling, as in the 
usual TLL theory. The key ingredient in the following 
is to consider quadratic ferromagnetic rather than lin- 
ear antiferromagnetic spin waves* 1 ^ The spin excitations 
are thus treated exactly as a Heisenberg ferromagnet re- 
siding on a fluctuating lattice, with the total Hamilto- 
nian decomposed into charge (H c ) and spin (H s ) parts, 
H = H c + H s , characterized by independent, sufficiently 
disparate propagating velocities, v c and v s , respectively. 
Although this type of separation of spin and charge exci- 
tations in ID Hubbard models, in the limit of very strong 
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repulsive interactions, is supported by exact Bethe ansatz 
calculations, 1 - further numerical work on more general 
ID and quasi- ID systems would be desirable to establish 
whether this property is more generic to strongly- coupled 
quasi-lD systems. Therefore, we proceed phenomenolog- 
ically, without providing rigorous microscopic justifica- 
tion, noting that other authors have developed similar 
approaches for ID antiferromagnetic systems^ii 

We shall demonstrate several salient features of G(x, t) 
in the regimes of J > and J = 0. When J = 0, 
the asymptotic form of G{x, t) is characterized by strong 
Gaussian decay over a distance, while retaining anoma- 
lous power-law singularities with interaction-dependent 
exponents. Furthermore, the low-energy behavior of the 
tunneling density of states retains the same exponent as 
the spin incoherent regime^ l/4g—l, where g is the TLL 
interaction parameter corresponding to a spinless chain 
with the same particle density and interactions. This 
implies that the system is also characterized by non-TLL 
universality. For J > 0, we find a peculiar enhancement 
of the plasmon peak spectral singularity that is sensi- 
tive to the spin kinematics. The effective coupling of the 
spin and charge collective modes is quantified in terms 
of a kinematic mismatch of the two coupled propagating 
excitations, which strongly suppresses the plasmon peak. 

The effective Hamiltonian. Let us begin with a chain 
of spin-5 particles, average interparticle spacing a and 
strong short-range repulsion. For the entirety of this pa- 
per, we shall work within the real time and space do- 
mains, as it will facilitate the proper physical interpreta- 
tion of our results. The charge sector retains the typical 
TLL form, 



dx 
2^ 
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where the bosonic <p c and C fields satisfy the usual com- 
mutation relations [9 c {x),(j) c {x 1 )] = — (in /2)sgn(x — x'), 
and g characterizes the renormalized interparticle repul- 
sion strength. The spin sector experiences the underlying 
fluctuating Wigner lattice with spins that interact via 
Heisenberg ferromagnetic exchange, 7 H s = — J^i^i^i ' 
Si+i, where I labels lattice sites, which are essentially 
distinguishable due to strong repulsion. The exchange in- 
teraction arises due to the remaining wave-function over- 
lap between the neighboring particles! 7 - To account for 
the lattice zero-point motion and vibrations, in the fol- 
lowing, the classical labels I will be replaced with the 
bosonic displacement operators. Furthermore, to be spe- 
cific, we set the nearest-neighbor exchange Ji = hJ /2ir 2 S 
to be constant, where J = irv s /a defines a characteris- 
tic spin velocity v s , and S is the particle spin in units 
of h. We note, however, that the exact form of the spin 
Hamiltonian should not be important, as long as we are 
dealing with a ferromagnetic ground state whose low- 
energy excitations are gapless magnons with quadratic 
dispersions. (It should also be noted that, because of 
the quadratic dispersion, there is no single spin velocity 
associated with the ferromagnetic magnons, unlike with 



the soundlike modes of the TLL.) Therefore, the long- 
wavelength free spin-wave propagator (on a stationary 
lattice) takes on the usual form: 



G, 



(S+(N,t)S-(0,0)) = Sy/{2ir/iJty^n- 



(2) 



Single-particle Green's function. In order to calculate 
the full Green's function, which reflects both spin and 
charge propagation, we have to take into account the 
fluctuations of our lattice. We shall make use of the first- 
quantized path-integral formalism to trace over the most 
relevant spin and charge degrees of freedom.— The num- 
ber of sites N(x,t) the spins have to propagate through 
is essentially the number of particles between (0, 0) and 
(x,t). N(x,t) fluctuates about a mean value fix, as fol- 
lows: N(x, t) =nx+ [8 c (x, t) - 6> c (0, 0)] /it. The central 
maneuver en route to computing a particular correlation 
function in a given regime is to account for the explicit 
dependence of N(x,t) on the bosonic charge degrees of 
freedom in the expectation brackets. 

We are interested in computing the single-particle 
Green's function G(x,t) = (ipi(x, t)ip\(0, 0))f t , for suffi- 
ciently large x and t, in the spin-up ferromagnetic ground 
state. (For spins larger than 1/2, by "spin down" we 
mean a spin S — 1 state.) This correlator character- 
izes the amplitude for an injected spin-down particle to 
propagate from to x in time t in a polarized spin-up 
environment. In order to properly account for charge 
correlations, we will employ the standard bosonization 
dictionary. The spinless bosonic field creation and anni- 
hilation operators are given by ip^ B (x,t) ~ e-^*'*) and 
ipB{x,t) ~ e i( t>e{x,t) ^ wnere the dual field 9 C is related to 
the particle density fluctuations at long wavelengths by 
p(x) = d x 6 c (x)/n. It is also necessary to account for the 
chirality of the excitations in terms of independent right- 
moving and left-moving parts, which can be conveniently 
labeled by the light-cone variables xl,r = x ± v c t. We 
shall ignore the time ordering associated with the factor 
[0 c (x), <j) c (x')], as it simply shifts the particle number by 
a constant factor. 

J = bosons. We now focus on the J = limit in 
a spin-up polarized ground state, where the spin-wave 
excitations are suppressed and only plasmons propagate, 
which was first discussed in Ref. 13. Since the particles 
cannot flip their spins, the injected spin-down particle be- 
comes effectively distinguishable and has to be removed 
in order to contribute to the propagator. This means 
that we have to force N(x, t) to vanish)^ 



Next, we expand the S function to obtain 



G(x,t)~ I —e 



(3) 



(4) 



where we have made use of the definitions <d(x,t) = 
9 c (x,t) - C (O,O), $(x,t) = Mx,t) - C (O,O), and 
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Furthermore, the Gaussian averages satisfy 



the identity (e / 



AA 2 }/2 



Therefore, after complet- 



ing the squares, performing the integral over A, and em- 
ploying the standard expressions for the bosonic cor- 
relators: (6 2 ) = g]n(xRX L ), (* 2 ) = {l/Ag)\n(x R x L ), 
($6) = (1/2) \ti{xr/xl), we finally obtain at long times 



G(z,i)~ — e 
u(t) 



XrX L 



(5) 



where u(t) = (a/ir)yj2g \n(v c t/a) is the root-mean- 
square particle displacement during a time t due to quan- 
tum fluctuations 32. One can easily notice that, because 
of the strong Gaussian decay on the scale of u, there is 
no appreciable contribution to the propagator close to 
the light cones x l , xr rs at large times (taking into ac- 
count a short-distance Luttinger-liquid cutoff ~ a). See 
the J > case below, however. 

J = fermions. We must account for the Pauli exclu- 
sion principle by the negative sign acquired by the many- 
body wave function when two fermions are exchanged. 
Consequently, the factor (— 1)^ must be inserted in the 
expectation brackets (J3]), but because the 5 function en- 
forces a vanishing N, for which this factor is inconsequen- 
tial, we obtain the same Green's functions ((5]). 

Clearly, the propagator ([5|) is dominated by a strong 
spatial Gaussian decay, and by inspection one can infer 
the asymptotic behavior of the local tunneling density of 
states v(e), recovering the familiar power-law singular- 
ity near the Fermi level that is reminiscent of the spin- 
incoherent regime^ v{e) ~ |e| 1 ^ 4ff_1 /-\/ln(l/ |e|), where 
the energy e is measured with respect to the Fermi level. 

J > bosons (S — 1). A ferromagnetic spin- wave 
texture is included by tracing over a product of the spin 
propagator @ with the standard TLL plasmon propa- 
gator. Operationally, a coupling of the two excitations 
is established by replacing the S function of Eq. with 
the spin propagator l|2]). in terms of the bosonized parti- 
cle number operator N: 



G(x,t) 



e ii3 ^ B (z,t)^(0,0) 



(6) 



Evidently, the iV 2 term gives rise to quadratic operators 
O 2 in the exponent that are hard to evaluate in gen- 
eral. Nevertheless, as a first step in understanding the 
basic physics, we shall employ the leading correction to 
the mean-field approximation, by retaining contributions 
linear in the &(x,t) fluctuations. This approximation 
should be valid at large x (see below for the exact con- 
dition), so that the fluctuations SN N. However, we 
concede that it is possible that the quadratic terms may 
provide meaningful corrections, and, in particular, would 
be required for a proper analytic continuation to imagi- 



nary time. In our approximation 
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governed by the exponents 
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2Jt 2 V Jt 
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J > fermions (S — 1/2). For the fermionic case, we 
follow similar steps, after inserting the additional factor 
(-1)^ = Re e mN = {e l7rN + e" i7rAr )/2 into the expecta- 
tion brackets. Let us denote G±(x,t) as corresponding 
to the e ±ivN parts, with G{x,t) = G+(x,t) +G-(x,t): 



G+(x,t) 



e ik F x+i F Jt ' 




(9) 



Adding the other G—(x, f) component yields 



G(x,t) ■ 
where 

Sr± 
Sl± 



e 2.it 



X R X, 



+ 



X R X, 



(10) 



— ± 1 

8g 



kFX 3/ k F x 

2Jt 2 V Jt 

kFX j/ k F x 

2Jt 2 I Jt 



(11) 



Notice that G-(x,t) = G+(—x,t), so that G(x,t) = 
G(—x,t), as required by the inversion symmetry. 

When J — > 0, Eq. ([6|) reduces to Eq. (|3|) and we do not 
have spin-wave propagation. In this regime, however, 
particle number N is forced to vanish and our mean-field 
treatment assuming SN <C N, which is required for the 
validity of Eqs. (|T))- (fTTj) . fails. It is important to derive 
the exact conditions for which our mean-field approxima- 
tion of the propagator dominates over the N ~ contri- 
bution (J5]). By comparing Eq. ([7]) with Eq. (O, we can 
see that when t < 1/J = l/kpv s = t*, the cx x 2 con- 
tribution to the exponents ([5]) actually leads to a more 
severe suppression of the Green's function than the Gaus- 
sian decay of Eq. ([5]) . [See the prefactor governed by the 
exponent Sl of Eq. (|12p below, which can be rewritten 
as the Gaussian exponential decay e~( 9lnt t 2 ^ x l Vst " 1 . The 
case of fermions is completely analogous.] Our mean-field 
analysis for a finite J will thus fail at t < f * , and the prop- 
agator will be dominated by the contribution (O decay- 
ing on the scale of the interparticle spacing a. For t ^> i*, 



4 




50 100 

Jt 

FIG. 1: Single particle Green's function G(x,t), Eq. (0, valid 
for Jt 3> 1. Here, we plot the absolute value (inserting a 
short-distance cutoff a: xr : l — > xr,l ± ia), with white color 
corresponding to zero. The dashed trace is the slice of the 
density plot at kpx = 50, while the solid trace also includes 
the oscillations of Re G(x, t) due to the spin propagator. Most 
visible features are due to the broad Gaussian peak near x — 
of width ~ v 3 t in the x direction and the plasmon peak 
near x = v c t [governed by the exponent ~ (g/2)(v c /v s ) 2 , see 
Eq. (j 1 2 p ] . Because the relative weight of the plasmon feature 
is suppressed with respect to the x = correlator by the 
factor ~ e -(9int/2)(v c /v s f ^ it wm Qnly be appreciable if the 

ratio Vc/vs is moderate. In the plot, we chose v c /v s — 3 and 
g — 1/2 (corresponding to a strong short-range repulsion). 
Note that the exponents (O approach l/8g at t 2> t s — x/v s . 

Eqs. (|7|)- (|lip are valid, corresponding to a Gaussian de- 
cay over a longer distance ~ v s t, where the spin propaga- 
tion is carried predominantly by spin waves rather than 
charge fluctuations. Zvonarev et al^ recently discussed 
the two-particle Green's function describing spin-density 
propagation, capturing the crossover from the trapped 
spin regime at J — > to propagating spin waves at a fi- 
nite J. They found that the crossover takes place when 
t ~ 1/J, irrespective of x, in agreement with our results. 

Discussion. It is clear that the asymptotic forms of 
G(x, t) discussed here do not fit within the conventional 
TLL form. When J = 0, there is a strong Gaussian decay 
of G(x,t) over relatively short distances [see Eqs. ©]r^ 
on which is superimposed the power-law decay reminis- 
cent of the spin-incoherent physics. 6 At a finite J and 
t > („ the Gaussian decay is replaced by spin oscilla- 
tions and some remnant plasmon features [see Eqs. l[7]). 

It is important that we attempt to clarify the physical 
meaning of the power-law exponents ([8]) associated with 
the plasmon features, for a finite J, assuming i » i,. 



(The case for fermions is very similar.) Let us fix the 
spatial coordinate x and vary time t > from about t* 
though and above the plasmon peak near xr ~ (which 
requires sufficiently large distances x ^> v c i* ^> a). The 
propagation of plasmon and spin- wave excitations are as- 
sociated with the characteristic time scales t c = x/v c 
and t s — x/v s , respectively, and for the consistency of 
our formulation, we require t s 3> t c . Retaining only the 
right-moving plasmon peak, we can approximate Eq. (JT]) 
as (omitting the obvious free spin propagator) 

G ~ (t - t c )- 5 * ~ e-^^\t - t c )-^ 2 , (12) 

assuming t s /t c 3> 1 and3I g ~ 1 (focusing on the range of 
t from about t* to some multiple of t c ). Green's function 
(Tj"2"]) has two salient features: Gaussian decay in space on 
the scale of v s t and a remnant power-law plasmon peak 
at t — > t c governed by a large exponent oc (t s /t c ) 2 . Let 
us try to understand the physical meaning of the ratio 
t s /t that is key to both features. First, we notice that 
the free spin- wave propagator, Eq. @, enters as a pref- 
actor in Eq. (|10p. When a particle is removed at the 
position x and time t, the phase of the spin propagator 
changes by ~ N/Jt — nt s /t. The ratio t s /t > 1 
thus corresponds to Aip > tt, which suggests that the 
fluctuations of N associated with the removal of a single 
particle lead to a destructive interference of the spin com- 
ponent of the full propagator, suppressing the Green's 
function. When t <C t s , the associated exponent scales 
as (t s /t) 2 , and when t ~ £*, this suppression becomes 
more severe than that of Eq. ([5]) , where spins are carried 
by the lattice fluctuations rather than than the spin flips. 
This reflects an amusing interplay of the spin and charge 
kinematics, which has no precise analogies in the con- 
ventional TLL theory. At much longer times t ^> t s , and 
fixed x, the asymptotic power-law exponent approaches 
8r + 3l — * 1 /4<7 and the zero-bias spectral anomaly is 
governed by the exponent l/4</— 1/2, similar to the spin- 
incoherent case^i A summary of this behavior is schemat- 
ically illustrated in Fig. [T] We finally note that the sup- 
pression of the free spin-wave correlator ([2]) by a factor 
~ e -9 ln *(W*) 2 reported in Ref. [13, at t > t*, appears 
to be also related to the "smearing" of the oscillating 
spin-wave propagator by the charge fluctuations, which 
is responsible for the suppression of the single-particle 
propagator (fl"2"|) . 

We are grateful to the authors of Ref. [l3| for sharing 
their results prior to publication as well as for useful dis- 
cussions. 
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